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In multitype lattice gas models with hard-core interaction of Widom-Rowlinson 
type, there is a competition between the entropy due to the large number of types, 
and the positional energy and geometry resulting from the exclusion rule and the 
activity of particles. We investigate this phenomenon in four different models on the 
square lattice: the multitype Widom-Rowlinson model with diamond-shaped resp. 
square-shaped exclusion between unlike particles, a Widom-Rowlinson model with 
additional molecular exclusion, and a continuous-spin Widom-Rowlinson model. In 
each case we show that this competition leads to a first-order phase transition at 
some critical value of the activity, but the number and character of phases depend 
on the geometry of the model. Our technique is based on reflection positivity and 
the chessboard estimate. 
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1 Introduction 



Although the most familiar examples of phase transitions in lattice models originate 
from a degeneracy of ground states and therefore occur at low temperatures, this is 
not the only situation in which phase transitions can occur. Another possible source 
of criticality is a conflict of energy and entropy. This was noticed first by Dobrushin 
and Shlosman || in the case of an asymmetric double-well potential with two (sharp 
resp. mild) local minima separated by a barrier. They found that for some specific 
temperature energy and entropy attain a balance leading to the coexistence of highl- 
and low-temperature phases corresponding to the two wells; cf. also Section 19.3.1 of ||. 
Later on, Kotecky and Shlosman [[K]] observed that such a first-order phase transition 
can occur even in the absence of an energy barrier, provided there is an "explosion" 
of entropy. They demonstrated this in particular on the prototypical case of the q- 
state Potts model on Z rf for large q, showing that for a critical temperature there exist 
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q distinct ordered low-temperature phases as well as one disordered high-temperature 
phase; see also Section 19.3.2 of ||. 

This paper has the objective of studying entropy-driven first-order phase transitions 
of similar kind in multitype lattice gas models with type-dependent hard core interac- 
tion. In such models the crucial parameter is the activity instead of temperature, and 
the entropy-energy conflict turns into a competition between the entropy of particle 
types and the positional energy and geometry resulting from the exclusion rule and the 
activity of particles. One is asking for a critical activity with coexistence of low-density 
and high-density phases. 

The basic example of this kind is the multicomponent Widom-Rowlinson lattice 



gas model investigated first by Runnels and Lebowitz [15] in 1974 and studied later 



(theoretically and numerically) by Lebowitz et al. [12], cf. also [|13|]. If the number q of 
types is large enough (the numerical estimates give q > 7), there exist three different 
regimes: besides the low-density uniqueness regime and a high-density regime with 
q "demixed" phases for z > z c (q), there exists an intermediate domain of activities 
z o(q) < z < z c (q) with two "crystal" (or "staggered") phases with an occupation pattern 
of chessboard type, and the phase transition at z c (q) is of first order. The transition 
between staggered and demixed phases at z c (q) is again entropy-driven: in the staggered 
phases the type entropy wins, with the effect of an entropic repulsion of positions forcing 
the particles onto a sublattice, whereas in the demixed phases the particles gain energy 
and positional freedom but loose their type entropy. The same kind of phenomenon 
has also been discovered for a class of spin systems with annealed dilution (including 
diluted Potts and plane rotor models) @, f|. 

The aim of the present paper is to analyze the interplay of type entropy and the 
geometry induced by the lattice and the exclusion rule. While we stick to the integer 
lattice Zi d (and for simplicity in fact to the case d = 2), we vary the exclusion rule 
in order to gain some insight into the geometric effects involved. We investigate and 
compare four different models: 

1. the standard multitype Widom-Rowlinson model; 

2. a multitype Widom-Rowlinson model with nearest-neighbor and next-nearest 
neighbor exclusion between particles of different type; 

3. a multitype Widom-Rowlinson model with additional type-independent hard-core 
interaction between next-nearest particles; 

4. a ferrofluid model of oriented particles with exclusion between neighboring but 
not sufficiently aligned particles. (Similar continuous-spin counterparts of models 
2 and 3 will also be considered.) 

We show that in each of these examples an entropy-driven first-order phase transition 
occurs, but the number and specific characteristics of coexisting phases are different 
in all cases. Our technique is similar to that used in ||, ^ ||, [l(|] and Chapters 18/19 
of (8|, and is based on the (trivial) reflection positivity in lines through lattice sites 
and the resulting chessboard estimate p|. While in model 1 this is only an alternative 
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(and perhaps more elementary) approach to the results obtained in [12] by means of 
Pirogov-Sinai theory, the very same argument works also in the other models with only 
slight modifications. 

This paper is organized as follows. In Section ^ we introduce the four models and 
present our results. The proofs follow in Sections || to |6[ The general scheme is explained 
in detail for model 1, the standard Widom-Rowlinson lattice model. In the other cases 
we only indicate the necessary changes. 

Acknowledgement. H.O.G. gratefully acknowledges warm hospitality of the Centre de 
Physique Theorique in Marseille-Luminy, and V.Z. of the Mathematical Institute of the Uni- 
versity of Munich. This work was supported by the Deutsche Forschungsgemeinschaft, SPP 
1033. 



2 Models and results 

2.1 The multitype Widom-Rowlinson lattice gas 

This model describes a system of particles of q different types ('colors') which are allowed 
to sit on the sites of the square lattice Z 2 . (For simplicity we stick to the two-dimensional 
case; an extension to higher dimensions is straightforward, cf. Chapter 18 of Q or 
Jt|.) At each lattice site we have a random variable cr, taking values in the set E = 
{0, 1, . . . , q}. The equality Oi = means that site i is empty, and <7j = a G {1, . . . , q} 
says that i is occupied by a particle of color a. Particles of different color interact 
by a hard-core repulsion: they are not allowed to sit next to each other. There is no 
interaction between particles of the same color. This means that the formal Hamiltonian 
has the form 

H(a)=J2u(* i ,<r j ), (1) 

where the sum extends over all nearest-neighbor pairs (ij) C Z 2 of lattice sites (i.e., 
\i — j\ = 1), and the potential U is given by 

U(a l ,a J ) = l C ° K0**i**J*0> (2) 
I otherwise. 

This model is a lattice analog of the continuum two-species model of Widom and Rowl- 
inson J16|| , and was first introduced by Lebowitz and Gallavotti (O] (for q = 2) and 
Runnels and Lebowitz |15[ (for general q). 

Since U is either or oo, the temperature does not play any role, and the only 
energetic parameter is the activity z > which governs the overall particle density; we 
assume that the activity does not depend on the particle color. Accordingly, the Gibbs 
distribution in a finite region A C Z 2 with boundary condition rj in A c = Z 2 \ A is given 
by 

= h^ v off A} (Z^)- 1 z N ^ exp [ - £ U(a u a,)} , (3) 

where N\(a) = \{i £ A : ex.; ^ 0}| is the number of particles in A, and Z^ q q is a 
normalizing constant. 
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Alternatively, we may think of as obtained by conditioning a Bernoulli measure 
on the set of admissible configurations. Let 

fl = {(r£ E z : V(ij) <Ji<Jj = or o~i = aj} 

be the set of all admissible configurations on Z 2 . Given any such configuration ae!! 
and any subset A of Z 2 , we write cja for the restriction of a to A. We also write CIati 
for the set of all admissible configurations a G E^ in A which are compatible with some 
7] G fi, in the sense that the composed configuration ar/\c belongs to Q. In particular, we 
write £l\ = ^a,o f° r the set of all admissible configurations in A, which are compatible 
with the empty configuration outside A. It is then easy to see that 

Va% = *a'(- l n A,.j) > 

where 7r^ 9 = ® igA n-' q is the A-product of the measures TT-' q = (jq^-j, ■ ■■ , j^) 
on E. 

Given the Gibbs distributions /i^'l, we define the associated class Q(z, q) of (infinite 
volume) Gibbs measures on fi in the usual way ||. Our main result below shows that 
for large q there exist two different activity regimes in which Q(z,q) contains several 
phases of quite different behavior. These regimes meet at a critical activity z c (q) and 
produce a first-order phase transition. 

The different phases admit a geometric description in percolation terms. Let Z 2 be 
equipped with the usual graph structure (obtained by drawing edges between sites of 
Euclidean distance 1). Given any a G a subset S of Z 2 will be called an occupied 
cluster if S is a maximal connected subset of {i G Z 2 : Oi ^ 0}, and an occupied sea if, in 
addition, each finite subset A of Z 2 is surrounded by a circuit (i.e., closed lattice path) 
in S. In other words, an occupied sea is an infinite occupied cluster with interspersed 
finite 'islands'. If in fact o~i = a for all i G S we say S is an occupied sea of color a. 
We consider also the dual graph structure of Z 2 with so-called *edges between sites of 
distance 1 or \/2, and the associated concept of *connectedness. An even occupied *sea 
is a maximal *connected subset of {i = (11,12) G Z 2 : i\ + i2 even, cij 7^ 0} containing 
*circuits around arbitrary finite sets A. Likewise, an odd empty *sea is a maximal 
*connected subset of {i = (11,12) G Z 2 : i\ + 12 odd, o% = 0} surrounding any finite A. 

Theorem 2.1 If the number q of colors exceeds some qo, there exists an activity thresh- 
old z c (q) G }q/5,5q[ and numbers < e(q) < 1/3 with s(q) — > as q — > 00 such that 
the following hold: 

(i) For z > z c (q), there exist q distinct translation invariant 'colored' phases fi a G 
Q(z, q), a G {1, . . . , q}. Relative to fi a , there exists almost surely an occupied sea of color 
a containing any given site with probability at least 1 — s(q). 

(ii) For qo/q < z < z c (q), there exist two distinct 'staggered' phases ii cvcn ,/U odd G 
Q(z,q) invariant under even translations. Relative to /U cvcn , there exist almost surely 
both an even occupied *sea and an odd empty *sea, and any two adjacent sites belong to 
these *seas with probability at least 1 — e(q). In addition, all occupied clusters are finite 
almost surely, and their colors are independent and uniformly distributed conditionally 
on their position. /u odd is obtained from /i cvcn by a one-step translation. 
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(iii) At z = z c (q), a first- order phase transition occurs, in the sense that q+2 distinct 
phases fi cvcn , ^ odd , \i\ , ... ,fj,g E G(z c (q),q) coexist which enjoy the properties above. 

The preceding theorem can be summarized by the following phase diagram. 

q + 2 

phases: ^ 1 2 staggered ^ q colored 

z c (q) 

We continue with a series of comments. 

Remark 2.1 (1) The existence of staggered phases in an intermediate activity region 
was first observed by Runnels and Lebowitz [15|. As will become apparent later, this is 



a consequence of the fact that the lattice Z is bipartite and the interaction is nearest 



neighbor. According to Theorem 2.1, for large q the staggered regime extends up to 
the fully ordered regime, and the transition from the staggered regime to the ordered 
regime at z c (q) is of first order. This result (which disproves a conjecture in fly!) has 



already been obtained before by Lebowitz, Mazel, Nielaba and Samaj [12|. While their 
argument relies on Pirogov-Sinai theory (which even gives the asymptotics of z c (q)), 
we offer here a different proof based on reflection positivity which is quite elementary 
and can easily be adapted to our other models (including a continuous-spin variant of 
the present model). 

(2) There are two kinds of ordering to be distinguished: positional order and color-order. 
The colored (or 'demixed') phases fix, . . . ,fj, q show color-order but no positional order. 
(The impression of positional order is a delusion coming from the lattice regularity.) 
Their high density takes advantage of the chemical energy of particles (i.e., of the activity 
z.) On the other hand, the staggered (or 'crystal') phases /i cvcn , /U odd exhibit positional 
order but color-disorder. Positional and color-disorder occurs in the uniqueness regime 
at sufficiently low activities. 

(3) The first-order transition at z c (q) manifests itself thermodynamically by a jump of 
the particle density as a function of the activity. In fact, z c (q) can be characterized as 
the unique value where the density jumps over the level 2/3, cf. Lemma |3.6| . 

(4) For small z there exists only one Gibbs measure in Q(z,q). For example, using 
disagreement percolation one easily finds that this is the case when qz < p c /(l ~ Pc), 
where p c is the Bernoulli site percolation threshold for Z 2 ; see [jl], ||] for more details. We 
do not know whether the uniqueness regime extends right up to the staggered regime. 
As will be explained in the next comment, this question is related to the behavior of 
the hard-core lattice gas. 

(5) As was already noticed in |l5|, |l2) , the occupation structure of the large-g Widom- 
Rowlinson model at activity z is approximately described by the hard-core lattice gas 
with activity £ = qz. This becomes evident from the following argument (which is more 
explicit than those in 15, 12]). Consider the Gibbs distribution n\. pcT of our model 



in a rectangular box A with periodic boundary condition. (This boundary condition 
is natural, since later on we will only look at phases appearing in the extreme decom- 
position of infinite volume limits of /UA, pcr ; we could also use empty or monochromatic 
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boundary conditions instead.) Let ipA ;Pei be the image of yUA, pe r under the projection 
a — > n = (rti)igA = (l{ CTi ^o})ieA from E A to {0, 1} A mapping a configuration of colored 
particles onto the occupation pattern. tp\ jPCI is called the site-random-cluster distribu- 
tion, see Section 6.7 of Q. Its conditional probabilities are given by the formula 

V>A, Pcr (^ = iK\ W ) = qz * z qK{iyn) , 

where k(i, n) is the number of clusters of {j G A \ {i} : rij = 1} meeting a neighbor 
of i. Now, since K(i,n) = if and only if all neighbors of i are empty, this conditional 
probability tends to the one of the hard-core lattice model with activity £ when q — > oo 
and (, = qz stays fixed. Unfortunately, this result is not sufficient to conclude that in 
this limit the transition point from the unique to the staggered phase converges to that 
of the hard-core lattice gas, although this seems likely and is suggested by simulations 

0- 

(6) One may ask whether the monotonicity of the transition from the staggered to 
the ordered regime can be deduced from stochastic monotonicity properties of the site- 
random-cluster model, as is possible in the Potts model. Unlike in the standard (bond) 
random-cluster model, such a stochastic monotonicity is not available Q. To obtain 
the existence of a unique transition point z c (q) we will therefore use the convexity of 
the pressure, which implies that the particle density is an increasing function of z. 

(7) As is often the case in this kind of context, our bounds on go are n °t very useful. 
They only allow us to conclude that we can take e.g. qo = 2 • 10 85 . Note that for small 



q the ordered regime still exists [15], but for q = 2 there is no staggered regime but 
instead a direct second order transition from the gas phase to the ordered phase. For 
this and more information about the minimal q at which the staggered phase appears 

see mm. 

2.2 The square-shaped Widom-Rowlinson lattice gas 

The standard Widom-Rowlinson model considered above is defined by the exclusion 
rule that no two particles of different color may occupy adjacent sites. Equivalently, 
one may think of the particles as having the shape of the diamond {x £ R 2 : ||x||i < 1}, 
and diamonds of different color are required to be disjoint. 

In this section we want to study a variant with different geometry: we identify a 
particle at position i 6 Z 2 with the suitably colored square {x G R 2 : ||a; — *||oo < 1}; anci 
we stipulate that squares of different color are disjoint, while squares of the same color 
may overlap. Alternatively, this assumption amounts to replacing Z 2 by its matching 
dual (Z 2 )*, which is obtained from Z 2 by keeping all nearest-neighbor bonds and adding 
bonds between diagonal neighbors of Euclidean distance y/2; for lack of a generally 
accepted name we call this lattice the face-crossed square lattice. Accordingly, we say 
that two sites i,j G Z 2 are *adjacent if \i — j\ = 1 or \^2, and we write (ij)* for such 
pairs of sites. Saying that squares of different colors are disjoint is then equivalent to 
saying that particles of different color do not sit on *adjacent sites. 

The point of considering this model is that the square-shape of particles fits better 
with the geometry of Z 2 than the diamond-shape in the standard Widom-Rowlinson 
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model. As a consequence, the somehow artificial staggered phases disappear, and the 
model is closer to what one would expect to hold in the continuum Widom-Rowlinson 
model. 

We now turn to precise statements. The formal Hamiltonian of the square-shaped 
Widom-Rowlinson model is analogous to (||), 

H(a) = J2 Ufa,*,) , 

for 

with the same pair interaction U as in (0). The variables <r, still take values in the set 



E = {0, 1, . . . , q}. We then can introduce the same definitions as in Section 2.1, only 
replacing (ij) by (ij)* at all proper places. For simplicity, we refrain from adding *'s 
to the quantities so defined. In particular, we write again Q(z,q) for the set of Gibbs 
measures in the present square-shaped model. 

Given a configuration a G $7, we call a subset S of Z 2 an occupied * cluster if S is 
a maximal *connected subset of {i 6 Z 2 : q / 0). Such an S may be visualized as 
the connected component in R 2 of the union of all squares with centers in S. We say 
that two disjoint square-shaped particles are contiguous if they touch each other along 
at least one half of a side (i.e., if their centers have Euclidean distance 2 or y/E), and 
are separated by two empty sites. A set S of pairwise disjoint square particles will be 
called a contiguity sea if it is maximal connected relative to the contiguity relation, and 
the union of their closures surrounds each bounded set. 

Our result for this model is the following. 

Theorem 2.2 // q exceeds some sufficiently large qo, there exist a critical activity 
z c(q) £ ] <? 1 ^ 3 /3, 3 q 1 / 3 [ and numbers < e{q) < 1/3 with e(q) — > as q — > oo such 
that the following hold: 

(i) For z > z c {q), there exist q distinct translation invariant 'colored' phases \i a G 
Q(z, q), a G {1, ...,<?}. Relative to fj, a , there exists almost surely an occupied sea of color 
a containing any given site with probability at least 1 — s(q). 

(ii) For qo/q < z < z c (q), there exists a translation invariant disordered phase 
^di S G Q{ z iQ) such that with probability 1 all occupied *clusters are finite, independently 
and randomly colored, and surrounded by a contiguity sea. Moreover, [i dis (o~i ^ 0) < 
1/4 + e(q) for all i. 

(iii) At z = z c (q), a first-order phase transition occurs, in the sense that there exist 
q + 1 distinct phases /i dis , fi\, . . . , \x q G Q(z c (q), q) exhibiting the properties above. 

In short, we have the following phase diagram: 

9 + 1 

phases: ^ 1 disordered ^ q colored 

z c (q) 

Remark 2.2 (1) The first-order transition at z c {q) manifests itself thermodynamically 
by a jump of the particle density from a value close to 1/4 to a value close to 1. 
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(2) The behavior of the square-shaped Widom-Rowlinson model differs from that of 
the standard, diamond-shaped Widom-Rowlinson model in that there are no staggered 
phases but instead only one disordered phase showing not only color-disorder but also 
positional disorder. In fact, we expect that this disordered phase is the unique Gibbs 
measure for any z < z c (q), although this does not follow from our methods. (Just as 
in Remark |2.1| (5), one can see that for large q the square-shaped Widom-Rowlinson 
model is related to the square-shaped hard-core lattice gas. It seems that the latter 
model does not exhibit a phase transition, but we are not aware of any proof.) 

(3) The first-order transition at z c (q) implies a percolation transition from an empty sea 
in the disordered phase to an occupied sea in the colored phases. In spite of the supposed 
uniqueness of the Gibbs measure in the whole range [0, z c (q)[, this interval contains a 
further percolation threshold, namely a critical value for the existence of a contiguity 
sea. Indeed, for sufficiently small z one can use disagreement percolation ||] to show 
the uniqueness of the Gibbs measure and the existence of a sea of empty plaquettes; 
the latter excludes the existence of a contiguity sea. However, it is far from obvious 
that a contiguity sea will set on at a well-defined activity. This is because neither 
the existence of a contiguity sea is an increasing event, nor the site-random cluster 
distributions (which can also be used in the present case) are stochastically monotone 
in z, cf. Remark |2.1| (6). We note that a similar percolation transition occurs also in 
the square-shaped hard-core lattice gas ||; this shows again that the latter describes 
the limiting behavior of our model in this regime. 

2.3 A Widom-Rowlinson model with molecular hard core 

Another way of changing the geometry of the Widom-Rowlinson model is to introduce 
a molecular (i.e., color-independent) hard-core interaction between particles. In this 
section we will discuss a model variant of this kind. 

As before, the underlying lattice is still the square lattice Z 2 , and the state space at 
each lattice site is the set E = {0, 1, . . . , q}. The formal Hamiltonian is of the form 

H(a)= £ *(«7 ij «7 i )+ E U(wi)\ ( 4 ) 
K-il=i \i-j\=V2 

here <J>, the nearest-neighbor molecular hard-core exclusion, is given by 

*(o-i,Oj-) -- 



cc if fjjCJj / , 
otherwise, 



and the next-nearest neighbor color repulsion U is still defined by (|2|). The main effect 
of the molecular hard core is a richer high-density phase diagram containing 2q phases 
with color-order and staggered positional order. The low-density regime is disordered 
both in the sense of color and position, as in the case of the square-shaped Widom- 
Rowlinson model. The transition between these regimes is still of first order, though 
the positional order of the high-density phases is an impediment for this to occur. We 
do not repeat the definitions of admissible configurations and of Gibbs measures, which 
are straightforward. 
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Theorem 2.3 // q > q for a suitable qg, there exist a threshold z c (q) G ]g/18, 18g[ 
and numbers < e(g) < 1/5 with e{q) — ► as q — > oo such that the following hold: 

(i) For z > z c (q) there exist 2q distinct colored and staggered phases /U ajCVCn , /U aj0dd G 

Q{z,q), a G {1, . . . ,q}, which are invariant under even translations. Relative to /i a , C vcn 

there exist almost surely both an even occupied *sea of color a and an odd empty *sea, 

and any two adjacent sites belong to these *seas with probability at least 1 — e(q). fi a ,odd 

is obtained from fi a ,c V cn by a one-step translation. 
1 /7 

(ii) For q ' jq < z < z c (q), there exists a translation invariant disordered phase 
Mdis £ G(z,q) such that with probability 1 all occupied *clusters are finite, independently 
colored with uniform distribution, and enclosed by a contiguity sea. Also, /i dis (o"j ^ 0) < 
l/4 + e(q) for all i. 

(iii) At z = z c (q), a first-order phase transition occurs, in the sense that there 
coexist 2q + l distinct phases // dis , /xi jCVC11 , m i0dd , ... , /ig, ovcn , /x 3i0dd G G(z c (q),q), with the 
properties above. 

We summarize this theorem by the following phase diagram: 



2q + l 

phases: ^ 1 disordered ^ 2q staggered colored 

z c (q) 

Remark 2.3 (1) The first-order transition at z c (q) manifests itself thermodynamically 
by a jump of the particle density from a value close to 1/4 to a value close to 1/2. 
(2) The Widom-Rowlinson model with molecular hard-core may be viewed as a com- 
bination of a lattice gas of hard diamonds and the square-shaped Widom-Rowlinson 
model. Its high-density regime inherits the staggered occupation pattern from the for- 
mer, and the color order from the latter. The effect of colors is still strong enough 
to produce a first-order transition, which is absent in the pure hard-diamonds model 



[ 1 2 1 . Just as in the case of Theorem 2.2, the low density regime is governed by the 
behavior of the square-shaped hard-core lattice gas. Indeed, it is not difficult to develop 
a random-cluster representation of the model and to show that its conditional proba- 
bilities converge to that of the square-shaped lattice gas when q — > oo but qz remains 
fixed; cf. Remark [O] (5). Mutatis mutandis, the comments in Remark |2.2| apply here 
as well. 

(3) One may ask what happens if we interchange the roles of $ and U in the Hamiltonian 
of Eq. (U), i.e., if there is a molecular hard core between diagonal neighbors of distance 
\/2, and a Widom-Rowlinson intercolor repulsion between nearest particles of distance 
1. In this case it is not hard to see that for any q > 2 and sufficiently large z there 
exist four different phases with positional order but color disorder. One of these phases 
(to be called the even vertical phase) has almost surely a sea of sites i = G Z 2 

which are occupied when i\ is even, and empty when i\ is odd. The other three phases 
are obtained by translation and/or interchange of coordinates. However, it seems that 
in this case the geometry of interaction does not exhibit the properties leading to a 
first-order transition, so that the transition to the low density regime is of second order. 
We will return to this point at the end of Section |^. 
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2.4 A continuous-spin Widom— Rowlinson model 

The multitype Widom-Rowlinson lattice model may be viewed as a diluted clock model 
for which each lattice site is either empty or occupied by a particle with an orientation in 
the discrete group of g'th roots of unity. This suggests considering the following plane- 
rotor model of oriented particles which may serve as a simple model of a ferrofluid or 
liquid crystal. 

Consider the state space E = {0} U S , equipped with the reference measure v = 
5q + A, where A is normalized Haar measure on the circle S 1 . As before, the equality 
<7i = means that site i is empty, while o~i = a £ S* 1 says that i is occupied by a 
particle with orientation a. The formal Hamiltonian is again given by (|l]), where the 
pair interaction U is now defined by 

jjt \ I oo if <Ji,(Jj £ S 1 , o, L ■ Oj > cos27ra , 
l ' J I otherwise 

for some angle < a < 1/4. This potential forces adjacent particles to have nearly 
the same orientation. The parameter a will play the same role as 1/q did before. The 
Gibbs distributions fi^t in a finite region A C Z 2 with boundary condition n and 
activity z > are defined by their densities with respect to the product measure 
which are again given by the right-hand side of equation (|3j). We write Q(z,a) for 
the associated set of Gibbs measures. Since U preserves the 0(2)-symmetry of particle 
orientations, the Mermin-Wagner-Dobrushin-Shlosman theorem (cf. Theorem (9.20) of 
H) implies that each such Gibbs measure is invariant under simultaneous rotations of 
particle orientations. 

Theorem 2.4 If a is less than some sufficiently small ao, there exist a critical activity 
z c (a) £ ] a~ 2 /18, 5a~ 2 [ and numbers < e(a) < 1/3 with e(a) — > as a — > such that 
the following hold: 

(i) For z > z c (a) there exists a dense 'ordered' phase fi old £ Q(z,a) exhibiting the 
translation invariance and O (2) -symmetry of the model. Relative to fi old , there exists 
almost surely an occupied sea containing any fixed site with probability at least 1 — e(a) 
(and on which the orientations of adjacent particles differ only by the angle 2ira). 

(ii) For Oq 1 < z < z c (a), there exist two distinct 'staggered' phases /i cvcn ,/U odd £ 
Q(z,a) which are invariant under particle rotations and even translations. Almost surely 
with respect to fj, even there exist both an even occupied *sea and an odd empty *sea, and 
any two adjacent sites belong to these *seas with probability at least l — e(q). In addition, 
all occupied clusters are almost surely finite, and conditionally on their position the 
distribution of orientations is invariant under simultaneous rotations of all spins in a 
single occupied cluster. fi odd is obtained from /U cvon by a one-step translation. 

(iii) At z = z c (a), a first-order phase transition occurs, in the sense that there exist 
three distinct phases /U even , /x odd , // old £ G(z c (a),q) with the properties above. 

The theorem above shows that the present model behaves similarly to the related finite- 
energy model considered in Presumably fi old is the unique Gibbs measure for z > 
z c (a), and there is a second-order transition from the staggered regime to the low- 
activity uniqueness regime. We thus have the following phase diagram. 
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phases: ^ 1 2 staggered ^ ordered 

z c (q) 

Remark 2.4 The model above is a continuous-spin counterpart of the standard 



Widom-Rowlinson model considered in Section |2.l| . It is rather straightforward to 
modify our techniques for investigating analogous continuous-spin variants of the square- 
shaped Widom-Rowlinson model and of the model with diagonal molecular hard core. 
In the first case, we obtain a phase diagram of the form 

2 

phases: 1 disordered ordered 



' z 



z c {q) 
and in the second case we find 



3 

phases: 1 disordered 2 staggered ordered 



' z 



z c (q) 
The details are left to the reader. 



3 Proof of Theorem [Z7[] 

The proof of all four theorems follows the general scheme described in Chapters 18 and 
19 of Georgii J8|], which is similar in spirit to that of Dobrushin and Shlosman [pi] and 



Kotecky and Shlosman [10]. This scheme consists of two parts: a model-specific contour 
estimate implying percolation of "good plaquettes" , and a general part deducing from 
this percolation the first-order transition and the properties of phases. We describe the 
general part first and defer the contour estimate to a second subsection. Many of the 
details presented here for the Widom-Rowlinson model carry over to the other models, 



so that for the proofs of Theorems 2.2 to 2.4 we only need to indicate the necessary 



changes. We note that our arguments can easily be extended to the higher dimensional 
lattices Z d using either the ideas of Chapter 18 of H or those of . 



3.1 Competition of staggered and ordered plaquettes 

We consider the standard plaquette C = {0, l} 2 in Z 2 as well as its translates C + i, 
i £ Z 2 . Two plaquettes C + i and C + j will be called adjacent if \i — j\ = 1, i.e., if C+i 
and C + j share a side. We are interested in plaquettes with a specified configuration 
pattern. Each such pattern will be specified by a subset F of the set of admissible 
configurations in C. For any such F we define a random set V(F) as follows. Let r\ 
and T2 be the reflections of C in the vertical resp. horizontal line in the middle of C, 
and r % = r^r 1 ^ the reflection associated to i = (ii,«2) £ Z 2 . We then let 

V(F) :a ^{ieZ 2 : r l a c+t G F} (5) 
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be the mapping associating with each a £ Q, the set of plaquettes on which a shows 
the pattern specified by F. (The reflections r % need to be introduced for reasons of 
consistency: they guarantee that two adjacent plaquettes may both belong to V(F) 
even when F is not reflection invariant, as e.g. the sets G ovcn and G odd below.) 

We are interested in the case when F is one of the following sets of 'good' con- 
figurations on C. These sets are distinguished according to their occupation pattern. 
Describing a configuration on C by a 2 x 2 matrix in the obvious way, we define 

• G stag = G evcn U G odd = {Q : 1 < a, b < q} U : 1 < a, b < q}, the set of all 
staggered configurations with 'diagonal occupations'. 

• G old = Ui<a<g G a = Ui< a <<j{(a a)}> ^ ne se ^ °^ an f un y ordered configurations with 
four particles of the same color. 

• G = G stag U G ord , the set of all good configurations. 

Our first objective is to establish percolation of good plaquettes, i.e., of plaquettes in 
which the configuration is good; the other plaquettes will be called bad. We want to 
establish this kind of percolation for suitable Gibbs measures uniformly in the activity 
z (provided z is not too small). A suitable class of Gibbs measures is that obtained by 
infinite-volume limits with periodic boundary conditions. 
For any integer L > 1 we consider the rectangular box 

A L = {-12L + 1,...,12L} x {-14L + 1,...,UL} (6) 

in Z 2 of size v(L) = 24 L x 28 L. (The reason for this particular choice will become 



clear in the proofs of Lemmas 3/? and [3^.) We write ^ z £ q per for the Gibbs distribution 
in Al with parameters z, q and periodic boundary condition, and Q pci (z,q) for the set 
of all limiting measures of fJ> z £ 9 peT as L — > oo (relative to the weak topology of measures). 
The basic result is the following contour estimate which shows that bad plaquettes have 
only a small chance to occur. 

Proposition 3.1 For any 5 > there exists a number q$ G N such that 

H(A n V(G) = 0) < <5 |A| (7) 
whenever q > go, zq > qo, \i S G pci (z, q), and A C Z 2 is finite. 

In the above, {A n V(G) = 0} is a short-hand for the event consisting of all a for which 
all plaquettes C + i, i S A, are bad; similar abbreviations will also be used below. 

The proof of the proposition takes advantage of reflection positivity and the chess- 
board estimate, cf. Corollary (17.17) of Q, and is deferred to the next section. We 
mention here only that qo is chosen so large that 

5{q) = q -im + g -i/i2 + ? -i/4 + ? -l/2 < s (8) 

when q > qo- It will be essential in the following that the contour estimate is uniform 
for z > 
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As an immediate consequence of the contour estimate we obtain the existence of a 
sea of good plaquettes. We will say that a set of plaquettes forms a sea if the set of 
their left lower corners is connected and surrounds each finite set. It is then evident 
that the existence of a sea of completely occupied plaquettes implies the existence of an 
occupied sea; likewise, the existence of a sea of plaquettes which are occupied on their 
even points implies the existence of an even occupied *sea. In this way, the concept of a 
sea of plaquettes is general enough to include all concepts of seas introduced in Section 
^. Specifically, for any F C Qc we define S(F) as the largest sea in V(F) whenever 
V{F) contains a sea, and let S(F) = otherwise. 

For z > qo/q we write O poi (z,q) for the set of all accumulation points (in the weak 
topology) of measures /x n G G pci (z n ,q) with z n — > z, z n > qo/q. The graph of the 
correspondence z — > peT (z, q) is closed; this will be needed in the proof of property 
(A2) below. 

Proposition 3.2 For any e > there exists a number q$ G N such that 

/x(0G 5(G)) >l-£ 
whenever fi G G P er{z,q), q > qo and z > qo/q. 

Proof. Note first that the contour estimate (M) involves only local events and therefore 
extends immediately to all \i G Q pci (z,q). The statement then follows directly from 



Proposition 3.1 together with Lemmas (18.14) and (18.16) of Q. The number 5 has to 
be chosen so small that 45(1 — 5<5) -2 < e . □ 

What is the advantage of having a sea of good plaquettes? The key property is that the 
sets G stag and G ord have disjoint side-projections. That is, writing b = {(0,0), (1,0)} for 
the two points on the bottom side of C we have 

u G G stag , a' G G ord dj) ^ <7j, , 

and similarly for the other sides of C. As a consequence, if two adjacent plaquettes 
are good then they are both of the same type, either staggered or ordered. Therefore 
each sea of good plaquettes is either a sea of staggered plaquettes, or a sea of ordered 
plaquettes. Hence 

{S(G) + 0} = {S(G st J + 0} U {S(G OId ) + 0} , 

and the two sets on the right-hand side are disjoint. Moreover, the sets G cvcn and 
G odd also have disjoint side-projections, and so do the sets G a , 1 < a < q. Therefore, 
the event {5(G stag ) ^ 0} splits into the two disjoint subevents {S'(G even ) ^ 0} and 
{S(G odd ) ^ 0}, and {S(G OTd ) / 0} splits off into the disjoint subevents {S(G a ) / 0}, 
1 < a < q. In other words, each sea of good plaquettes has a characteristic occupation 
pattern or color corresponding to a particular phase, and we only need to identify the 
activity regimes for which the different phases do occur. 

To this end we fix any e > 0. We will need later that e < 1/6. As in the proof of 
Proposition 3.2, we choose some < 5 < 1/25 such that 45(1 - 55)~ 2 < e , and we let 
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qo be so large that condition (||) holds for all q > q^. For such q and q we consider the 
two activity domains 

Artag = [z > qo/q ■ M° ^ V(G Btag )) > fi(0 G V(G old )) for some n G £ pcr (z,g)} 

and 

Ard = [z > qo/q : M° G V(Gw)) > /u(0 G F(G stag )) for some /u G Q pci (z, <?)} . 

Our next result shows that these sets describe the regimes in which staggered resp. 
ordered phases exist. The mean particle density g(fi) of a measure /u which is periodic 
under translations with period 2 is defined by 

g(fi) = fi(N c )/\C\ ; (9) 

recall that Nq is the number of particles in C. 

Proposition 3.3 (a) For each z G ^4 stag there exist two 'staggered' Gibbs measures 
Aievenj/^odd ^ G ( z , q) invariant under even translations ofZi 2 and permutations of particle 
colors. n cvcn -almost surely we have S(G even ) ^ 0, and all occupied clusters are finite and 
have independently distributed random colors. In addition, ^ cvcn (0 G S(G cvcn )) > l — 2e, 
and in particular g(^ cvcn ) <\ + £. ^ Q dd has the analogous properties. 

(b) For each z G ^4 ord there exist q 'colored' translation invariant Gibbs measures 
Ha G G(z,q), a G {1, ■■■,(?} • Each fi a satisfies fi a (S(G a ) ^ 0) = 1, ^ a (0 G S(G a )) > 
1 — 2e, and in particular has mean particle density g{[i a ) > 1 — 2e. 

Proof, (a) Let z G ^4 stag be given and fi G Gp CI (z,q) be such that ^(0 G V(G sta _ g )) > 
fi(0 G V(G ord )). Then /x(0 G V(G old )) < 1/2 and therefore 



//(O G S(G st j) = M (0g5(G), 0^(G o 



1 1 

^ 1 - £ -2 = 2- £> °- 

But G stag splits into the two parts G cvcn and G odd which are related to each other by the 
reflection in the line {xi = 1/2}, and \i is invariant under this reflection. Hence 

p = /i(o G S(G cvcn )) = fi(0 G S(G odd )) > i Q - > . 

We can therefore define the conditional probabilities /i cvcn = n(- \S(G cvcn ) ^ 0) and 
A'odd = /■*( ' I'S'C^odd) 7^ 0)- Since the events in the conditions are tail measurable, these 
measures belong to Q(z,q). It is clear that these conditional probabilities inherit all 
common invariance properties of \x and the conditioning events. Moreover, we find 



(S(G evea ) ^ 0) = ~/x(S(G stag )^ 

< ~/x(OGS(G stag ))^ M (o£S(G) 

< P +|, 



14 



and therefore 

/w(o G S(G avcn )) > P > 1 - 2e . 
\ ' p + e/2 

In particular, it follows that 

£>(/4vcn) < ^ Meven(0 G y(G cvcn )) + /i cvcn (0 F(G cvcn )) < - + E . 

Finally, we show that /i cvcn -almost surely all occupied clusters are finite, and their colors 
are conditionally independent and uniformly distributed when all particle positions are 
fixed. Indeed, since jU cvon (5(G evcn ) ^ 0) = 1 there exists // cvcn - almost surely an odd 
empty *sea. This means that any box A is almost surely surrounded by an empty 
*circuit. On the one hand, this shows that all occupied clusters must be finite almost 
surely. On the other hand, for any r/ > we can find a box A' D A containing an 
empty *circuit around A with probability at least 1 — rj. Let T be the largest set with 
A C r C A' such that there are no particles on its outer boundary dT; if no such set 
exists we set T = 0. The events {r = A} then depend only on the configuration in Z 2 \A. 
By the strong Markov property of // ovcn , we conclude that on {r ^ 0} the distribution 
of colors of the occupied clusters meeting A is governed by the Gibbs distribution in T 
with empty boundary condition. The symmetry properties of the latter thus imply that 
these colors are conditionally independent and uniformly distributed. Letting r] — * 
and A | Z 2 we find that this statement holds in fact for all occupied clusters. 

By construction, fi odd is obtained from /i even by a one-step translation, and thus has 
the analogous properties. 

(b) The proof of this part is quite similar. Pick any z G A OTd and \x G pei (z, q) such 
that n(0 G V(G old )\0 G V(G)) > 1/2. Since /i is invariant under permutations of colors 
it then follows in the same way that 

P^(0G5(G fl )) >IQ_ e ) > , 

so that we can define the conditional probabilities [i a = fi(-\S(G a ) ^ 0) G Q(z,q), 
a £ {1, ... ,q}. Also, 

n(s(G a ) ^ 0) < i /i(0 G S(G ord )) + %(0 S(G)) <P+ £ -> 

whence /i a (0 G S{G a )) > p/(p + e/q) > 1 - 2s and g(ji a ) > fi a (0 G V{G a )) > 1 - 2e. □ 

According to the preceding proposition, Theorem will be proved once we have 
shown that there exists a critical activity z c (q) G ] q/5, 5q [ such that ^4 stag = [qo/q, z c (q)] 
and A OTd = [z c (q), oof. To this end we will establish the following items: 

(Al) A stag U A old = [qo/q, oof . 

(A2) A atag and A OTd are closed. 

(A3) A old n [qo/q, q/5] = 0. 

(A4) ^ stag n [5g,oo[ = . 
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(A5) |^4 stag n A old \ < 1 . 

Statement (Al) follows trivially from the definitions of A stag and vl rd- Assertion (A2) 
is also obvious because these definitions involve only local events, and the graph of the 
correspondence z — ► Q pei (z, q) is closed by definition. 

Property (A3) corresponds to the discovery of Runnels and Lebowitz [15] that stag- 
gered phases do exist in a nontrivial activity regime, and follows directly from the next 
result. 

Lemma 3.4 For z < q/5 and \x G Q pci (z, q) we have 

Jo G V(G ord ) G V(G)) < 1/2 . 



Proof. Consider the Gibbs distribution /j>l in the box Al with periodic boundary 
condition, and let 

Gord,L = [o- G O ijPer : a- C {i) G G ord for all i G A L | ; (10) 

here we write f2z, iPcr for the set of admissible configurations in the torus A^ (including 
nearest-neighbor bonds between the left and the right sides as well as between the top 
and bottom sides of A^), and C(i) for the image C + i mod Al of C under the periodic 
shift of Al by i. (As G ord is reflection-symmetric, we can omit the reflections r l which 
appear in (||).) The chessboard estimate (cf. Corollary (17.17) of M) then implies that 

We compare the latter probability with that of the event 

G eveni L = |er G ^L.pcr : r % a C (i) G G cvcn for all i G A L \ . (11) 

This gives 

j£lJG ord ,L) < VlU G °^)/JlU G ^l) = X« L \ z -«M q -v(L)/2 

because G old< L contains only the q distinct close packed monochromatic configurations, 
while for a G G eveUj L the v(L)/2 particles can have independent colors. Taking the 
v(L)'th root and letting L — > oo we find for \i G G pai (z,q) 

M (0 G V{G OTd j) < (z/qfl 2 < 5- 1 / 2 < (1 - S)/2 . 

The last inequality comes from the choice of 5. Since ^(0 G V(G)) > 1 — 5 by Proposition 
3.1, the lemma follows. □ 

Assertion (A4) corresponds to the well-known fact that q ordered phases exist when 
the activity is large. For q = 2 this was already shown by Lebowitz and Gallavotti p[ , 



and for arbitrary q by Runnels and Lebowitz [15|. This is again a simple application of 
the chessboard estimate. 

Lemma 3.5 For z > 5q and fj, G G peT (z,q) we have 

fi(o G ^(G s ta g )|o G V(G)j <l/2. 
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Proof. Let G ord ,L be as in (|i"ol) , and define G stag) £ analogously. By the chessboard 
estimate we find 

— ( / i L,pcr(^'stag,L) / Ml ,p 

(G 

ord,_L/ / / 

< 2 lML) z l/2 q l/2 z -l q -l/v(L) 

because G atag) L = G CVC[1) l U C? oddj £ contains 2 q v ( L )l 2 distinct configurations of particle 
density 1/2. We can now complete the argument as in the preceding proof. □ 

For the proof of (A5) we will use a thermodynamic argument, namely the convexity 
of the pressure as a function of log z. For any translation invariant probability measure 
fx on Q we consider the entropy per volume 

s(fi)= Km IAI-^Gua). 

|A|^oo 



Here we write //a for the restriction of fi to Q 



A, 



is the entropy of /^a, and the notation |A| — ^ cxo means that A runs through a specified 



increasing sequence of square boxes; for the existence of s(/x) we refer to [g, |14 |. 
We define the thermodynamic pressure by 

P(log z) = max g(/_t) log z + s(fi) ; (12) 

L 

the maximum extends over all translation invariant probability measures \x on S7, and 
q(/j,) = /x(oo 7^ 0) is the associated mean particle density, cf. @. (Since $7 is defined 
as the set of all admissible configurations, the hard-core intercolor repulsion is taken 
into account automatically.) By definition, P is a convex function of logz, and the 
variational principle (see Theorems 4.2 and 3.12 of [p[| ) asserts that the maximum 



in (12) is attained precisely on Ge(z,q), the set of all translation invariaion invariant 
elements of Q(z,q). By standard arguments (cf. Remark (16.6) and Corollary (16.15) 
of Q) it follows that P is strictly convex, and 

PLQogz) < q(ji) < P+(logz) for all p G Ge(z,q) ; (13) 

here we write P'_ and P' + for the left-hand resp. right-hand derivative of P. By strict 
convexity, P'_ and P' + are strictly increasing and almost everywhere identical. Assertion 
(A5) thus follows from the lemma below. 

Lemma 3.6 For each z G A stag n A OTd we have P'_ (log z) < 2/3 < P+(logz). 



Proof. This has already been shown essentially in Proposition |3.3j . Pick any z G A stag n 
A ord , and let \x G Q pci (z,q) be as in the proof of Proposition |3.3| (a). Consider the 
conditional probability /i stag = /i( • |5(G stag ) ^ 0) = \ /U oven + \ /x odd . By the arguments 
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there, /i stag is well-defined, belongs to Ge(z,q), and satisfies g(fJ- stag ) < 1/2 + e < 2/3. 
On the other hand, the measures \x a constructed in Proposition |3.3| (b) also belong to 
Ge(z, q) and satisfy g(fJ, a ) > 1 — 2e > 2/3. The lemma thus follows from (|l3|). □ 



We can now complete the proof of Theorem 2.1. Properties (Al) to (A4) together 
imply that A sta _ g n A old ^ 0. This is because the interval [qo/q,oo[ is connected and 
therefore cannot be the union of two disjoint non-empty closed sets. Combining this with 
(A5) we find that A stag n A old consists of a unique value z c (q). In particular, A old cannot 
contain any value z < z c (q) because the infimum of such z's would belong to A sta _ g n A old ; 
likewise, A sta _ g does not contain any value z > z c {q). Hence A sta _ g = [qo/q, z c (q)] and 



A old = [z c (q), oo[, and Theorem 2.1 follows from Proposition 3.3. 



3.2 Contour estimates 



In this subsection we will prove Proposition Consider the set Q,c of all admissible 
configurations in C, and the set B = Qc \G of all bad configurations in C. We split B 
into the following subsets which are distinguished by their occupation pattern: 

• Bq = {(„ jj)}, the singleton consisting of the empty configuration in C. 

• B 1 = { (° °) , (° °) , (° p) , (q °) : 1 < a < qr}, the set of all configurations with a single 
particle in C. 

• £ 2 = O, O, O : 1 - a - ?}' the set of admissible configurations for 
which one side of C is occupied, and the other side is empty. 

• B 3 = (°°) : 1 < a < g}, the set of all admissible configurations 
with three particles in C. 

We then clearly have B = Ufc=o Bk ■ The four different kinds of "badness" of a plaquette 
will be treated separately in the three lemmas below. We start with the most interesting 
case of plaquettes with three particles. 
For any L > 1 and k G {0, . . . , 3} let 

Bk,L = {c G O ijPeI : a c{i) G S fc for all i G A L } , 

where C{i) is as in (ll0|). Consider the quantities p z k ,q L = ^, z £ q peT (Bk^) 1 / v ^ and p z k ,q = 
limsuPL-^ooPjfe',*- 

Lemma 3.7 < g -1 / 56 /or all z > and q G N. 

Proof. Fix any integer L > 1 and consider the set -E^l of configurations a in A^ having 
a single empty site in each plaquette. We claim that \B^l\ < q 2 14i+2 . First of all, for 
each a G B^l the occupied sites in A^ form a connected set, so that all particles have 
the same color. Thus there are only q possible colorings, and we only need to count 
the possible occupation patterns for a G B%^. It is easy to see that the plaquettes 
C(i) with a(i) =0 form a partition of A/,. For each such partition, the plaquettes are 
either arranged in rows or in columns. In the first case, each row is determined by 
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its parity (even or odd), namely the parity of i\ for each C{%) in this row; likewise, in 
the second case each column is determined by its parity. We can therefore count all 
such partitions as follows. There are 4 possibilities of choosing the plaquette containing 
the origin. If this plaquette is fixed, there are no more than 2 14L ~ 1 possibilities of 
arranging all plaquettes in rows and choosing the parity of each row. Similarly, there 
are at most 2 12L_1 possibilities of arranging the plaquettes in columns. The number of 
such partitions is therefore no larger than 4(2 14L ~ 1 + 2 12L ~ 1 ), and the claim follows. 

To estimate fj, z £ q per (B 3 ^) we will rearrange the positions of all particles so that many 
different colors become possible. More precisely, we divide into (3L)(4L) rectangular 
cells A(j) of size 8x7. Let Ao(i) be the rectangular cell of size 7x6 situated in the 
left lower corner of A(j), and consider the set 

F 3 ,l = We n LyPCI : a + on A (j)> ° = on A(j) \ A (j) for all j} . 

Since |A(j) \ Ao(j)| = 8-7 — 7-6 = |A(j)|/4 for all j, each a £ F^l has particle 
number 3t>(L)/4, just as the configurations in -63,1,. As the colors of the particles in 
the blocks Ao(j') can be chosen independently, we have \F 3: l\ = q 12L2 = q v ( L ^ 56 . (The 
above construction, together with a similar construction in the proof of the next lemma, 
explains our choice of the rectangle A^.) Now we can write 

,z,q 1 r> \ ^ ^L 9 P cS B 3,l) \B 3tL \ 14L+2 l-„(L)/56 

The proof is completed by taking the v(L)'th root and letting L — ► 00. □ 

Next we estimate the probability of plaquettes with two adjacent particles at one 
side of C . 

Lemma 3.8 p z 2 ,q < q~ 1/12 for all z > and qeN. 

Proof. Fix any L > 1, and let a G B 2 ,l- Then the particles are either arranged in 
alternating occupied and empty rows, or in alternating occupied and empty columns. 
The colors in all rows resp. columns can be chosen independently of each other. Hence 
\B 2 , l\ = 2(g ,14i + q 12L ) < Aq 14L . Moreover, each a G B 2: l has particle number v(L)/2. 
As in the last proof, we construct a set F 2> l of configurations with the same particle 
number but larger color entropy as follows. 

We partition A^ into (8L)(7L) rectangular cells A(j) of size 3x4, and let Ao(j) be 
the rectangular cell of size 2 x 3 in the left lower corner of A(j). We then define 

F 2 ,l = {a£ Q LtPCI : a + on A (j), o = on A(j) \ A (j) for all j} . 

Since | A(j) \ Ao(j')| = 3 • 4 — 2 • 3 = | A(j')|/2 for all j, each a G F 2 ^ has particle number 
v{L)/2. As the particle colors in the blocks Ao(j) can be chosen independently, we have 
1-^3, L I = g 56 ^ 2 = q v ( L y 12 . As in the last proof, we thus find 

I R A ^ VlIcSB^l) \B 2 ,l\ . n UL-v(L)/12 
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Taking the v (L)'th root and letting L — > oo we obtain the result. □ 

Finally we consider the probability of 'diluted' plaquettes with a single or no particle. 

Lemma 3.9 p^ q < (zq)^ 1 / 2 and p\' q < (zq)' 1 / 4 for all z > 0, q G N. 
Proof. We consider first the case of no particle. For each L > 1 we can write 

z,q lTD _ N / VL q pe X B 0,L) _ 1 



/4? pcr (5 ,L) ^ 



/x^ pei (G evett)i ) ' 

where G evcnj L is defined by (|Tl|). The identity follows from the facts that -Bo,L contains 
only the empty configuration, whereas each configuration in G evcnj L consists of v(L)/2 
particles with arbitrary colors. The first result is thus obvious. 

Turning to the case of a single particle per plaquette, we note that each a G B\l 
consists of v(L)/4 particles with arbitrary colors, and there are no more than 2 14i+2 
distinct occupation patterns for these particles; the latter follows as in the proof of 



Lemma 3.7 (by interchanging empty and occupied sites). Hence 



» L , pel {Xl,L) _ ^ {Qc _ l) - z v(L)/2 qV (L)/2 > 

and the second result follows by taking the u(L)'th root and letting L — > oo. □ 



Proof of Proposition 3.1. Let \i G Q pcr {z, q) and a finite A C Z 2 be given. Then we can 
write 

/x(A n V(G) = 0) = ^ : G for all i G A) 

7 :A^{0,...,3} 

< limsu P/ i L,pcr(°' : a C(i) G s 7 (i) for all i G A) 

7 :A^{0,...,3} L ^°° 
7:A^{0,...,3} L ^°° iGA 



< (l>r) 

\ I A / 



|A| 



•fc=0 

In the third step we have used the chessboard estimate, see Corollary (17.17) of |§. 



Inserting the estimates of Lemmas 3/7 , 3^ and 3J3 and choosing qo as in (||) we get the 
result. □ 



4 Proof of Theorem 12.2 



Here we indicate how the proof of Theorem 2.1 can be adapted to obtain Theorem |2.2| . 
First of all, the different geometry of the present model leads to a new classification 
of good and bad plaquettes: the ordered configurations in G old are still good, but the 
(former good) configurations in G stag are now bad and will be denoted by -B stag , while 
the configurations in B\ are now good, and we set G dis = B\. 
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We first need an analog of the contour estimate, Proposition [O]. Remarkably, the 
estimates of Lemmas 3/7 and 3^ carry over without any change. To deal with £> stag 
we can proceed exactly as in Lemma 3.8, noting that each configuration in -B stagi L is 



monochromatic, so that |-B stag ,L| = 2q. This shows that also Pst2 g < q 1 ^ 12 - Finally for 
Bq we compare the set -Bq,l with 



Fi, L = {ae fl L , pcr : Oi £ iff i G 2 Z 2 } ; 



(14) 



this gives p z ' q < (zq) 1//4 . The counterpart of Proposition |3j] thus holds as soon as 



is so large that 



-1/56 



+ 2, 



-1/12 



-1/4 



< 5. 



With the contour estimate in hand we can then proceed as in Section |3.1| . Proposition 
3.3 carries over verbatim; the only difference is that G stag is replaced by G dis (which is 
not divided into two parts with disjoint side-projections), and g(^ diB ) < \(l — 2e) + 2e = 
z + ^f- By the latter estimate, the assumption e < 1/6 is slightly stronger than necessary 



for adapting Lemma 3.6 to the present case, but we stick to it for simplicity. 



The counterparts of Lemmas 3.4 and 3.5 are obtained as follows. On the one hand, 
we have the estimate 

/4',UG OTd ,L) < ^JG otd!L )/rtlJF l>L ) < Z^q z -^)/ V ,(L)/4 

showing that 

fi(0 G V(G ovd )) < (z 3 A) 1/4 < 3~ 3/4 < (1 - 8)12 
when [i G G pci (z, q) and z < g 1 / 3 /3. On the other hand, as in Lemma |3.9| we find 

2 l4L+2 z v(L)/4 q v(L)/4 



and therefore 



ML). 



fi(0 G V(G dis )) < (g/z 3 ) 1/4 < 3" 3 / 4 < (1 - 5)/2 
when [i G G pC r(z,q) and z > 2>q l l^. With these ingredients it is now straightforward to 



complete the proof of Theorem 2.2 along the lines of Section 3.1 



5 Proof of Theorem 1273 



Here we consider the Widom-Rowlinson model with molecular hard-core exclusion. We 
look again at good configurations in plaquettes. The set £Iq of admissible configurations 
in C splits into the good sets 

G ord = G 

even 

l<a<q 1<<J<<? 

of ordered staggered configurations, the good set 

G dis =£i^ax^o,a 

of disordered configurations, and the only bad set Bq consisting of the empty config- 
uration. The main technical problem which is new in the present model is that the 
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sets G ord and G dis fail to have disjoint side-projections (although this is the case for the 
sets G^cven and G a ,odd)- We therefore cannot simply consider sets of good plaquettes, 
but need to consider the sets of "good plaquettes with neighbors in the same phase". 
Accordingly, we introduce the random sets 

V(G old ) = {ieV(G old ):i + (l,0),i + {0,l)eV(G old )}, 
V(G dis ) = {iey(G dis ):i + (l,0),i + (0,l)ey(G dis )}, 

and V(G) = V(G old ) U V^(G dis ). By definition, a sea in V(G) then contains either a 
sea in V(G old ) or a sea in V(G dis ). To establish the existence of such a sea we use the 
following contour estimate. 

Proposition 5.1 For any S > there exists a number qo G N such that 

n(AnV(G) =0) < <5 |A| 

whenever q > qo, zq > q^ 7 , [i G G pcr (z, q), and A C Z 2 is finite. 

Proof. Let us start by introducing some notations. We consider the sublattices 

£i,even = {i = (ii,«2) G Z 2 : h is even} , £i j0dd = Z 2 \ £i, ovon , 

and their rotation images C2, avcn and £2,odd which are similarly defined. We also intro- 
duce the horizontal double-plaquette 

D 1 = CU(C + (1,0)) ={0,1,2} x {0,1} 

and the event 

Ei = jcr G E Dl : a c G G md , (?c+(i,o) G ^dis, or vice versaj 

that the two sub-plaquettes of D\ are good but of different type. E\ thus consists of 
the configurations of the form (q ^ q) with 1 < a < q, and their reflection images. In 
the same way, we define the vertical double-plaquette D 2 = C U (C + (0, 1)) and the 
associated event E 2 . With these notations we have 

Z 2 \ V(G) c\JW k , 
k=i 

where the random subsets W k of Z 2 are given by 

W 1 = V(B ), W 2 = V(B ) - (1,0) , W 3 = V(B )-(0,1), 
W4 = {i G £i,even : o Dl+i G E ± } , W 5 = {i G £i j0dd : a Dl+i G Ei} , 
W 6 = {i G £ 2 ,ovcn : o- D2+i G E 2 } , W 7 = {i G £ 2 ,odd : o"D 2 +i G #2} • 

(The sets W k are not necessarily disjoint.) So, for each /x G (? pe r(z, (?) we can write 

A i(AnV r (G) = 0)< mma(A k cW k ), (15) 

A 1 U...UA 7 =A- fe - 7 
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where the sum extends over all disjoint partitions of A. We estimate now each term. 



Consider first the case k = 1. Just as in Lemma p.9| we obtain from the chessboard 
estimate 

z,q i d \i/v(L) 

where Fi^ is given by fll^). The same estimate holds in the cases k = 2, 3 because these 
merely correspond to a translation. 

Next we turn to the case k = 4. Let L be so large that D A4. Using reflection 
positivity in the lines through the sites of £i )CVCn , we conclude from the chessboard 
estimate that 

^ pcr (A 4 C W.f'^ < tilJ E ^) 2ML) 

for the event 

Ei,l = {cr G fiL, P cr : VD^i) G E\ for all i e A L n £i, cvcn } . 

In the above, D\{i) stands for the image D\+i mod of D\ under the periodic shift 
by i of the torus A/,. Each a 6 £i,z has the following structure: every fourth vertical 
line (with horizontal coordinate either or 2 modulo 4) is empty, and on each group 
of three vertical lines between these empty lines every second site is occupied, with 
the coordinates of occupied sites being either even-odd-even in these three lines, or 
odd-even-odd; see the figure below. 



o • o • o o 

o o • o o • 

o • o • o o 

o o • o o • 

o • o • o o 

o o • o o • 

o • o • o o 



Of course, the interaction implies that the color of particles is constant in each of these 
groups of three vertical lines. Consequently, each such a has particle number 3v(L)/8, 
and \E 1>L \ = 2(2q) 6L ; recall the definition (§) of A L . 



We now make a construction similar to that in Lemma 3.7. We divide Al into 12L 
rectangular cells A(j) of size 8x7. Let Aq(j) be the rectangular cell of size 7x6 
situated in the left lower corner of A(j), and consider the set 

Fl = {cr E ^L, pcr : <7j 7^ iff ii + Z2 is even and i £ Ao(j) for some j} . 

Since |A(j)\Ao(i)| = for all j, each a £ Fl has particle number 3v(L)/8, just 

as the configurations in Ei l. As the colors of the particles in the blocks Ao(j) can be 
chosen independently, we have = q 12Ij2 = q^( L )/ 56 . Hence 

,z,q s ^ ^L 9 paI (El,L) \Ei L \ 6L v ( L )/ 56 

^ JEl ' L) - » l Ufd ~ W - 2(2q) q 

and therefore, by taking the 2/-u(L)'th power and letting L —* 00, we obtain 

M (A 4 C ^ 4 ) 1/|A41 < q- 1/28 . 
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The same estimate holds in the cases k = 5, 6, 7, as these are obtained by a translation 
or interchange of coordinates. 

We now combine all previous estimates as follows. Let qo be so large that 
7 ((fa 1//28 ) 1 / 7 < S, and suppose that q > qo and zq > q^ 7 . Then /x(Afe C Wfc) < g "' A ' : ' / ' 28 
for all and thus, in view of (|l5|) and since |A&| > |A|/7 for at least one k, 



/.(A n v(G) = 0) < £ 1/28 ) |A|/r < ^ |A| - 

AiU...UA 7 =A 



The proof of the contour estimate is therefore complete. □ 



To prove Theorem 2.3 we can now proceed as in Section 3.1. Let S(G) be the largest 
sea in V(G) if the latter contains a sea, and S(G) = otherwise. It is then immediate 
that a counterpart of Proposition |3.2| holds, and the definition of V(G) implies that 



{5(G) + 0} = {S(G dis ) + 0} U {S(G md ) + 0} , 
where the two sets on the right-hand side are disjoint. Moreover, 

{S(G old ) + 0} C |J {S(G a , cvcn ) + 0} U {S(G a , odd ) / 0} . 

a=l 



By the argument of Proposition 3J3 we thus obtain the existence of 2q ordered phases (as 
described in Theorem |2.3| (i)) whenever z is such that /j(0 G V(G old )) > fi(0 G V(G diB )) 
for some \i G ^ pcr (-z, <?), and the existence of a disordered phase /j> dia whenever the reverse 
inequality holds for such a fi. We have g(fi diB ) < j + ^f and e(fJ. a , CV cn) = Q(^a, dd) > | — £• 
The topological argument of Section |3.1| together with obvious counterparts of Lemmas 



3.4 and [D] then show that both cases must occur simultaneously for some z = z c (q), 
and this z is unique by the convexity argument of Lemma |3.6| . (For the latter we need 
to assume that e < 1/10.) 

We conclude this section with a comment on the model with nearest-particle color 
repulsion and a molecular hard-core exclusion between next-nearest neighbors. 

Comment on Remark |2.3| (3). If the roles of $ and U are interchanged, the good 
ordered configurations in C are those with two particles of the same color on one side 
of C, and no particle on the opposite side; we call this set again G ord . For large z, 
one can easily establish a contour estimate implying the existence of a sea S(G OTd ), and 
thus by symmetry also the existence of the four phases mentioned in Remark (3). 
The disordered good plaquettes are again described by the set G dis . As in the case of 
the Hamiltonian (|j), the sets G ord and G dis have no disjoint side-projections. However, 
whereas in that case we were able to show an entropic disadvantage in having adjacent 
G ord - and G dis -plaquettes, this is not true in the present case. The configurations re- 
sulting from iterated reflections of a double plaquette of type (°°°) have the maximal 
entropy possible for this particle number. Therefore the system can freely combine or- 
dered and disordered plaquettes, and our argument for a first-order transition breaks 
down. So it seems likely that the transition from the ordered to the disordered phase is 
of second order. 
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6 Proof of Theorem 



The analysis of the plane-rotor Widom-Rowlinson model is very similar to that of the 
standard Widom-Rowlinson model; only a few modifications are necessary. We define 
again the set flc of admissible configurations in the plaquette C in the obvious way, 
introduce the sets G stag = G C¥C11 UG odd as in Section |3j] (replacing {1, . . . , q} by S 1 ), and 
set G ord = {a G Sic : Cj G 5* 1 for all i E C} and G = G ord U G stag . The main task is to 
obtain a counterpart of the contour estimate, Proposition 3.1. To this end we consider 



the same classes B^, k G {0, . . . , 3} of bad configurations as in Section (with the 
obvious modifications), and the sets B^^l and the associated quantities p z k ' Q . 

To deal with the case k = 3 we proceed as in Lemma |3.7|, arriving at the inequality 



z,a m s < / J L, P cr(- B 3,L) _ v A {B^ L ) 

Now, v^{B% l) < 2 14:L+2 (2a) 3v ( L ^ 4: ~ 1 ; the first factor estimates the number of possible 
occupation patterns, and the second term bounds the probability that the configuration 
is admissible (by keeping only the bonds in a tree spanning all occupied positions). On 
the other hand, v^{F^^l) > (a 7 ' e ~ 1 ) v ( L ^ 5G , as can be seen by letting the spins in each 
block Aq(A;) follow a "leader spin" up to the angle 2ixaj2. Hence 



< 2 14L+2 2 3u(£)/4 a v(L)/56-l 



and therefore p z 3 ' a < 2 3 / 4 a 1 / 56 . 

In the case k = 2 we proceed as in the proof of Lemma |3.8| . On the one hand, 

v K {B 2 ,l) < 2((2a) 24L - 1 ) 14L + 2((2a) 28L " 1 ) 12L < 4(2a)^ 2 - 14L 

since the spins are ordered in separate rows or columns, and 2a < 1. On the other 
hand, v^{F2,l) — (a 2 ' 3- 1 ) V ( L " 12 by the same argument as above. Hence 

(B 2 L ) < < 4 • 2^/ 2 a v ^' 12 (2ay UL 

and therefore p z ^ a < 2 1 / 2 a 1 / 12 . 
Finally, for fc = we obtain 



r&JBod < ^ JGcvcnL) - zV(L)/ 2 

and thus Pq 01 < z~ l l 2 . Likewise, in the case k = 1 we get as in Lemma |5! 



nT (Bil) 2 14L+2 ^( l )/ 4 



/*K.(<W) - 

and thereby p\' a < z _1//4 . Combining these estimates as in the proof of Proposition 3.1 
we arrive at the counterpart of ([?]) as soon as ao is so small that 



25 



and a < otQ, z > 1/a.Q. 



To complete the proof of Theorem |2.4| as in Section |3.l| we still need to adapt Lemmas 



3.4 and 3.5. Writing 



t* L>pe A^,L) _ ^J GcvcnL) ~ z v(L)/2 

we find that for z < a~ 2 /18 and [i G Q pai (z, a) 

fifo e y(G ord )) < z^ 2 2a < 2/y/TE < (1 - S)/2 . 



Likewise, since 



^L.porl^stag^J S z,a /(-i \ — 



we see that for z > 5 aT 2 and /i G Q pcT {z, a) 

ti(0 e V(G st j) < z- l ' 2 a- 1 < (1 - S)/2 . 



The remaining arguments of Section 3.1 can be taken over with no change to prove 



Theorem 2.4 
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